Tunable singlet-triplet splitting in a few-electron Si/SiGe quantum dot 
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We measure the excited-state spectrum of a Si/SiGe quantum dot as a function of in-plane mag- 
netic field, and we identify the spin of the lowest three eigenstates in an effective two-electron 
regime. The singlet-triplet splitting is an essential parameter describing spin qubits, and we extract 
this splitting from the data. We find it to be tunable by lateral displacement of the dot, which 
is realized by changing two gate voltages on opposite sides of the device. We present calculations 
showing the data are consistent with a spectrum in which the first excited state of the dot is a 
valley-orbit state. 



Silicon quantum dots are candidate hosts for semicon- 
ductor spin qubits, both because of long spin relaxation 
and coherence times for electrons in Si, and because of 
potential synergy with classical microelectronics. Long 
spin relaxation times have been demonstrated in Si quan- 
tum dots and donors [iHl], and measurements of ensem- 
bles of donor-bound spins by electron spin resonance have 
demonstrated T2 coherence times up to two seconds [5]. 
One of the key properties of silicon quantum dot spin 
qubits is the ability to tune in real-time tunnel rates and 
couplings between neighboring dots by controlling elec- 
trostatic gate voltages [6-8 . Tunable, gate-defined sili- 
con quantum dots often are designed to sit at the inter- 
face between pure Si and a barrier of either SiGe [H [10] 
or Si02 IIH13]. 

For initialization and readout of singlet-triplet (S-T) 
spin qubits, an essential parameter is the energy dif- 
ference between the singlet and triplet states of 
two electrons in one dot [HI [15]. The energy ^st is 
equal to the lowest single-particle excited state energy, 
less a correction arising from electron-electron interac- 
tions. In silicon nanostructures, which have states aris- 
ing from the two low-lying valleys in the Si conduction 
band, the sharpness and quality of the interface between 
Si and the SiGe or Si02 barrier material play an impor- 
tant role in determining this energy [16]. Experiments 
have shown that quantum-confined structures can have 
reasonably large valley splitting, ranging from 100 /ieV 
to 1 — 2 meV [TTHTOj . The existence of large valley split- 
tings in silicon quantum dots have led to large ^st and 
the observation of Pauli spin blockade [20l[2T]. However, 
systematic control of the valley splitting or, more directly, 
^ST has not been demonstrated in a silicon quantum dot. 

In this letter, we report a magnetospectroscopy study 
of a Si/SiGe double quantum dot with 2 and valence 
electrons on the left and right dots, respectively. We 
use a pulsed-gate voltage technique to measure the evo- 
lution of the ground and low-lying excited states of the 
left dot as a function of an in-plane magnetic field B. 
We extract the magnetic field 5st at which the ground 
state changes from singlet to triplet, corresponding to 
the Zeeman energy equaling the singlet-triplet splitting 
E'sT for the (2,0) charge configuration. We find that 
is tunable by lateral displacement of the quantum 




FIG. 1. (a) Scanning electron micrograph of a double dot 
identical to the one used in the experiment. The transition 
measured here is between the left dot (white circle) and the 
left reservoir. Charge sensing is performed by measuring the 
current /qpc through the QPC channel under a bias of 500/iV. 
(b) Stability diagram of the double dot with effective electron 
occupation numbers labeled. The white symbols between re- 
gions (1,0) and (2,0) correspond to the gate voltages for the 
data reported below in Fig. 3. The transition line at the 
bottom right of the plot is invisible, because the tunnel rate 
between the right reservoir and the right dot is very slow in 
that gate voltage regime. 

dot location, achieved by simultaneously tuning voltages 
applied to two gates on opposite sides of the dot. 5st 
evolves systematically as a function of the gate voltages, 
and we measure a fractional change in 5st of up to 19%. 
Changes in gate voltages can alter both the position and 
shape of the electron wavefunctions in quantum dots [22r - 
24 . Applying asymmetric changes to the voltages on ei- 
ther side of the quantum dot, as we do here, will change 
primarily the position of the quantum dot. We perform 
calculations showing that the fractional change in 5st 
we observe is consistent with valley-orbit mixing arising 
from a rough Si/SiGe interface, and that a change in po- 
sition alone is sufficient to account for the magnitude of 
the observed changes in 5st- 

A double quantum dot, shown in Fig. [ija), is fabri- 
cated as described in [2]. A quantum point contact is 
defined by gates R and Q and is used to perform charge 
sensing measurements. Gate L is connected to a pulse 
generator (Tektronix AFG 3252B), allowing the applica- 
tion of fast voltage pulses. The dc gate voltages are tuned 
so that the double dot is in the few-electron regime, as 
shown in Fig. [TJb). The change in background grayscale 
arises from changes in the QPC sensitivity caused by ca- 
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pacitive cross-talk in the device [25]. Measurements are 
performed in a dilution refrigerator at an electron tem- 
perature Te = 143 ±10 mK, determined as described in 
Ref. [2 . The electron occupation numbers are effective; 
we believe there are spin-zero closed shells of electrons in 
both the left and the right dots that do not participate 
in the physics discussed here. 

We determine the 2-electron singlet-triplet (S-T) split- 
ting by using the charge sensing quantum point contact 
to measure the gate voltage dependence of the transition 
to the 2-electron state as a function of B. Figs. |2|a) 
and (b) show the transconductance Gl = S^qpc/^^L as 
a function of measured with a lock-in amplifier using 
a 120 /iV ac voltage applied to gate L. The bright peak in 
the color plot corresponds to adding one electron to the 
left dot. The gate voltage of this transition first increases 
and then decreases as a function of B. 

The electrochemical potential /iN, and equivalently the 
gate voltage of transitions like those in Figs.[2|a) and (b), 
has a dependence on the in-plane magnetic field of the 
form djii^/dB = gfiB^Stot{N) [26 . Here g is the Lande 
^-factor, /iB is the Bohr magneton, and AStot{N) is the 
change in the z component of the total spin when the 
A^th electron is added to the dot. The electrochemical 
potential has a slope of +^/iB/2 when a spin- up electron 
is added (magnetic moment anti-parallel to 5), whereas 
the addition of a spin-down electron results in a slope 
of —g 11^/2 (magnetic moment parallel to B). The pos- 
itive slope in Fig. [2](a,b) at small B corresponds to the 
addition of a spin-up electron, forming a 2-electron spin- 
singlet ground state. At the value of B marked with the 
arrows in panels (a) and (b), the slope changes; for B 
larger than this value, the added electron is spin-down, 
and the ground state is the triplet T_. As indicated 
schematically in Fig. |2|c), the turning point of the slope 
corresponds to a magnetic field 5st at which the Zee- 
man energy of the state T_ is equal to Est{B = 0). The 
value of BsT is different in Figs.[2|a) and (b), indicating 
that Est is tunable with gate voltage. We discuss the 
physics of this tunability in detail below. 

Pulsed-gate spectroscopy methods [271 HI] allow us to 
confirm the state identification described above, while si- 
multaneously mapping out the excited-state energy spec- 
trum as a function of B. Square wave voltage pulses of 
peak-to-peak amplitude 3.6 mV and frequency 50kHz are 
applied to gate L, and the time- averaged value of Gl is 
recorded as a function of B, as shown in Fig.[2|d). Here, 
the bottom (top) line corresponds to the positive (neg- 
ative) edge of the pulse bringing the 2-electron ground 
state into resonance with the Fermi level of the lead. 
Both of these lines therefore reproduce the shape of the 
line in Fig.[2];b). 

The two middle lines in Fig. [2|d) meet at 5 = and 
correspond to the triplet states T_ and Tq, which are 
degenerate at this point, and as B increases, the lines 
split. The T_ line has a negative slope, corresponding 
to the addition of a spin-down electron, and this state 
becomes the ground state when B = Bst- The Tq line 
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FIG. 2. (a), (b) Ground state magnetospectroscopy for two 
different sets of gate voltages chosen so that the gate voltages 
for (b) favor a dot position farther to the right than those for 
(a) (see Fig. |3|. The positive slope for small B corresponds 
to a singlet ground state, and the negative slope at large B 
corresponds to a triplet T_ ground state. The turning point 
of the transition occurs when the Zeeman shift for the T_ is 
equal to the zero-field ^st- (c) Schematic diagram showing 
the transition as a function of B. (d) Excited state magne- 
tospectroscopy using pulsed-gate voltages for the dot position 
corresponding to (b). Transitions to the three lowest states, 
S, T_, and To are clearly visible in the figure, (e) Simulated 
excited-state magnetospectroscopy for the data in panel (d). 

has positive slope, corresponding to the addition of a 
spin-up electron. The triplet T+ state does not appear, 
as the addition of a single spin to the spin-down ground 
state cannot create a state with Sz = +1 . 

Figure 2(e) shows a theoretical simulation of the exper- 
iment of Fig.[2](d), performed using a coupled rate equa- 
tion model similar to that described in the supplemental 
material for Ref. [2 . The model includes thermal broad- 
ening but neglects energy-dependent tunneling. The S, 
To, and T_ loading and unloading rates and the temper- 
ature are determined by fitting the simulation to the data 
in Fig. ^d). We find the loading rates =45.1 kHz, 
= 216 kHz, and F^^ = 377 kHz, and the unloading 
rates Fg = 164 kHz, F^_ = 354 kHz, and F^^ = 183 
kHz. 

Using the method illustrated in Figs.[2|a) and (b), we 
measure the transition field Bst at four different gate 
voltage configurations, corresponding to the triangle, cir- 
cle, square, and diamond shown in Fig. [TJb). Along this 
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FIG. 3. BsT, the magnetic field at which the ground state 
shifts from singlet to triplet, for different sets of gate voltages 
Vl and Vr, corresponding to the symbols on the stability 
diagram in Fig. 1(b). Error bars are determined by the un- 
certainty in linear fits to lines like those in Fig. 2(a) and (b). 
The three gray lines correspond to three different sets of mi- 
croscopic parameters that can be used to fit the experimental 
results. 

line in the stability diagram, changes in and Vr tend 
to shift the dot physically from left to right as (Fl) 
is made more positive (negative). As shown in Fig. [sj we 
observe a systematic increase in Bst as we move from 
left to right in the stability diagram. Over this range, 
Bst increases from 1.68 ± 0.09 T to 2.02 ± 0.07 T, a total 
increase of about 19 %. 

The energy of singlet-triplet splitting can be expressed 
as Est = Ei — Eq Cqi — Cqo + ^st, where Eq and Ei 
denote the ground and first excited-state energies, Cqi 
and Coo are the Coulomb interaction energies of the two 
electrons in the singlet and triplet states, and Kst is the 
exchange energy [29 . A simple shift of the dot position 
is expected to have little effect on the last three terms, 
which correspond to interactions between the electrons. 
Similarly, the shape of the wavefunction envelope should 
change very little as a function of the dot displacement. 
Thus, it is important to check whether microscopic fea- 
tures of the quantum device can account for the changes 
in Est that we observe. 

The single-particle spacing AE = Ei—Eq has a contri- 
bution arising from the difference in valley components 
in the two lowest lying orbital states. The quantum well 
interface will have atomic steps and other sharp changes 
in potential that vary as a function of lateral position, 
and these variations can modify the coupling of the two 
z-valleys, contributing to a position dependence of the 
energy difference A£: [30l|3T]. 

To test whether a small atomic-scale variation can ac- 



count for the magnitude of the observed variations in 
Bst, we perform tight-binding calculations of the single 
particle energy levels of an electron confined near a sin- 
gle atomic step. The calculations use a two-dimensional 
tight-binding Hamiltonian similar to Refs. [32, 33 , in- 
cluding a parabolic lateral confinement potential. The 
fitting procedure varies the position of the atomic step, 
the parabolic confinement length scale, and the vertical 
electric field, enabling a calculation of the variation in 
as a function of gate voltage. To compare with the 
measured ^st, we also fit the sum of the Coulomb and 
exchange energies (Coi — Cqo + ^st), and the results 
are plotted in Fig. |3] The fitting is underconstrained, as 
there are many physical ways to produce similar valley 
splitting. To indicate the types of variations possible, 
three results are plotted in Fig. |3] as the solid, dashed, 
and dotted lines, and all three calculations can repro- 
duce the magnitude of the observed change in Bst- The 
the lowest excited-state can be classified as "orbital-like" 
when the calculated wavefunction contains a lateral node 
or "valley-like" when it does not, and both cases occur. 
For valley-like excitations, lateral translation of the dot 
with respect to a step results in a tunable valley splitting. 
For orbital-like excitations, strong valley-orbit coupling 
enables a tunable orbital energy splitting [34J. 

In conclusion, by performing magnetospectroscopy in 
a double quantum dot in Si/SiGe at the (1,0) to (2,0) 
charging transition, we measure the singlet-triplet split- 
ting for the (2, 0) charge configuration and find it to be 
tunable by displacement of the dot. Theoretical calcu- 
lations show that atomic-scale structure of the quantum 
well interface is sufficient to produce valley-orbit mixing 
large enough to account for the experimental observa- 
tions, even under the assumption that the dot shape is 
unchanged as a function of gate voltage. The change in 
Est is large enough to be easily measured and offers an 
additional control knob for manipulation and initializa- 
tion of singlet-triple quantum dot spin qubits. 

We acknowledge experimental assistance from B. Rose- 
meyer and D. Greenheck. This work was supported 
in part by ARO and LPS (W911NF-08-1-0482), NSF 
(DMR-0805045), and the United States Department of 
Defense. The views and conclusions contained in this 
document are those of the authors and should not be 
interpreted as representing the official policies, either ex- 
pressly or implied, of the US Government. This research 
utilized NSF-supported shared facilities at the University 
of Wisconsin-Madison. 



[1] M. Xiao, M. G. House, and H. W. Jiang, Phys. Rev. 
[ Lett. 104, 096801 (2010) 

[2] C. B. Simmons, J. R. Prance, B. J. Van Bael, T. S. Koh, 
Z. Shi, D. E. Savage, M. G. Lagally, R. Joynt, M. Friesen, 
S. N. Coppersmith, and M. A. Eriksson, Phys. Rev. Lett. 



106, 156804 (2011). 
[3] R. R. Hayes, A. A. Kiselev, M. G. Borselh, S. S. Bui, E. T. 
Croke, P. W. Deelman, B. M. Maune, I. Milosavljevic, J.- 
S. Moon, R. S. Ross, A. E. Schmitz, M. F. Gyure, and 
A. T. Hunter, "Lifetime measurements (Ti) of electron 



4 



spins in Si/SiGe quantum dots," ArXiv:0908.0173. 

[4] A. Morello, J. Pla, F. Zwanenburg, K. Chan, K. Tan, 
H. Huebl, M. Mottonen, C. Nugroho, C. Yang, J. van 
Donkelaar, A. Alves, D. Jamieson, C. Escott, L. Hollen- 
berg, R. Clark, and A. Dzurak, Nature 467, 687 (2010). 

[5] A. M. Tyryshkin, S. Tojo, J. J. L. Morton, H. Riemann, 
N. V. Abrosimov, R Becker, H.-J. Pohl, T. Schenkel, 
M. L. W. Thewalt, K. M. Itoh, and S. A. Lyon, "Electron 
spin coherence exceeding seconds in high purity silicon," 
ArXiv:1105.3772vl. 

[6] N. M. Zimmerman, B. J. Simonds, A. Fujiwara, Y. Ono, 
Y. Takahashi, and H. Inokawa, Appl. Phys. Lett. 90,1 
n~ 033507 (2007) 

[7J C. B. Simmons, M. Thalakulam, B. M. Rosemeyer, 
B. J. V. Bael, E. K. Sackmann, D. E. Savage, M. G. 
Lagally, R. Joynt, M. Friesen, S. N. Coppersmith, and 
M. A. Eriksson, Nano Lett. 9, 3234 (2009). 

[8] L. A. Tracy, E. P. Nordberg, R. W. Young, C. B. Pinilla, 
H. L. Stalford, G. A. T. Eyck, K. Eng, K. D. Childs, 
J. R. Wendt, R. K. Grubbs, J. Stevens, M. P. Lilly, M. A. 
Eriksson, and M. S. Carroll, Appl. Phys. Lett. 97, 192110 
(2010). 

[9] K. A. Slinker, K. L. M. Lewis, C. C. Haselby, S. Goswami, 
L. J. Klein, J. O. Chu, S. N. Coppersmith, R. Joynt, R. H. 

^Blick, M. Friesen, and M. A. Eriksson, New J. Phys. 7j 

24 6 (2005) 

[lOJ T. Berer, D. Pachinger, G. Pillwein, M. Miihlberger, 
H. Lichtenberger, G. Brunthaler, and F. Schaffler, AppL] 



r~Phys. Lett. 88, 162112 (2006) 

[Tl] S. J. Angus, A. J. Ferguson, A. S. Dzurak, and R. G. 
Clark, Nano Lett. 7, 2051 (2007) 

[12] E. P. Nordberg, H. L. Stalford, R. Young, G. A. T. Eyck, 
K. Eng, L. A. Tracy, K. D. Childs, J. R. Wendt, R. K. 
Grubbs, J. Stevens, M. P. Lilly, M. A. Eriksson, and 
M. S. Carroll, Appl. Phys. Lett. 95, 202102 (2009). 

[13] M. Xiao, M. G. House, and H. W. Jiang, Appl. Phys. 
Lett. 97, 032103 (2010)^ 

[14] J. Levy, Phys. Rev. L ett. 89, 147902 (2002)| 

[15] J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, 
A. Yacoby, M. D. Lukin, C. M. Marcus, M. P. Hanson, 
and A. C. Gossard, [Science 309, 2180 (2005) 

[16] M. Friesen, M. A. Eriksson, and S. N. Coppersmith, 
Appl. Phys. Lett. 89, 202106 (2006) 

[17] S. Goswami, K. A. Slinker, M. Friesen, L. M. McGuire, 
J. L. Truitt, C. Tahan, L. J. Klein, J. O. Chu, P. M. 
Mooney, D. W. van der Weide, R. Joynt, S. N. Copper- 
smith, and M. A. Eriksson, Nat. Phys. 3, 41 (2007), 

[18] M. G. Borselh, R. S. Ross, A. A. Kiselev, E. T. 



[19 

[20; 

[21 
"[22 

[23; 
[24 
[25; 



[26 
[27 



[28; 



[29; 



[30 
r — 

[31 

[32; 
[33; 

[34 



Croke, K. S. Holabird, P. W. Deelman, L. D. Warren, 
L Alvar ado- Rodriguez, L Milosavljevic, F. C. Ku, W. S. 
Wong, A. E. Schmitz, M. Sokolich, M. F. Gyure, and 
A. T. Hunter, Appl. Phys. Lett. 98, 123118 (2011). 
W. H. Lim, F. A. Zwanenberg, C. H. Yang, and A. S. 
Dzurak, "Spin filling of valley-orbit states in silicon quan- 
tum dot," ArXiv:1103.2895. 

N. Shaji, C. B. Simmons, M. Thalakulam, L. J. Klein, 
H. Qin, H. Luo, D. E. Savage, M. G. Lagally, A. J. Rim- 
berg, R. Joynt, M. Friesen, R. H. Blick, S . N. Copper- 
smith, and M. A. Eriksson, Nat. Phys. 4, 540 (2008) 
H. W. Liu, T. Fujisawa, Y. Ono, H. Inokawa, A. Fuji- 
wara, K. Takashina, and Y. Hirayama, Phys. Rev. B TT] 
_073310 (2008) 

J. Kyriakidis, M. Pioro-Ladriere, M. Ciorga, A. S. 
Sachrajda, and P. Hawrylak, Phys. Rev. B 66, 035320 
(2002). 

S. Amasha, K. MacLean, L P. Radu, D. M. Zumbuehl, 
M. A. Kastner, M. P. Hanson, and A. C. Gossard, Phys.j 
Rev. Lett. 100, 046803 (2008) 

M. Thalakulam, C. B. Simmons, B. M. Rosemeyer, D. E. 
Savage, M. G. Lagally, M. Friesen, S. N. Coppersmith, 
and M. A. Eriksson, Appl. Phys. Lett. 96, 183104 (2010). 
C. B. Simmons, M. Thalakulam, N. Shaji, L. J. Klein, 
H. Qin, R. H. Bhck, D. E. Savage, M. G. Lagally, S. N. 
Coppersmith, and M. A. Eriksson, Appl. Phys. Lett. 91, | 
_213103 (2007) 

Y. Hada and M. Eto, Phys. Rev. B 68, 155322 (2003). 
J. M. Elzerman, R. Hanson, L. H. W. van Beveren, 
L. M. K. Vandersypen, and L. P. Kouwenhoven, A ppl.] 
Phys. Lett. 84, 4617 (2004) 

M. Thalakulam, C. B. Simmons, B. J. V. Bael, B. M. 
Rosemeyer, D. E. Savage, M. G. Lagally, M. Friesen, S.N. 
Coppersmith, and M. A. Eriksson, Phys. Rev. B 84, 
045307 (2011). 

R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha, 
_and L. M. K. Vandersypen, Rev. Mod. Phys. 79, 1217| 
(2007) 



M. Friesen and S. N. Coppersmith, Phys. Rev. B 81, 
115324 (2010) 

D. Culcer, X. Hu, and S. Das Sarma, Phys. Rev. B 82, 
205315 (2010). 

T. B. Boykin, G. Klimeck, M. A. Eriksson, M. Friesen, 
S. N. Coppersmith, P. von Allmen, F. Oyafuso, and 
S. Lee, Appl. Ph ys. Lett. 84, 115 (200 4) 
A. L. Saraiva, B. Koiller, and M. Friesen, Phys. Rev. B 
82, 245314 (2010). 

J. K. Gamble, S. N. Coppersmith, and M. Friesen, To 
be submitted. 



